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1. INTRODUCTION 

Electricity plays an important role in all sides of human life. Therefore the increase of the world 
energy demand, due to the population growth, modern industrial society and the environment pollution, is 
moved the world towards renewable energy sources as the solution of these issues related to the energy 
demand, high fuel cost, and greenhouse problems, to enhance power quality issues and energy efficacy[1]. 

Out of all renewable energy units, wind and solar systems are considered as secured and reliable 
sources and being installed widely. These generation units have obtained popularity because of the 
environment-friendly characteristics and they are inexhaustible energy sources as well as fast development in 
the technologies [2]. The use of these sources recently attended a salient increase according to the countries' 
development and environmental pollutions. The generation powers from this renewable generation are 
intermittent that causes some mismatches between produced power and demanded one. This affects the 
microgrid operation and causes instability in the frequency and power of the hybrid energy system [3]. To 
tackle this problem, some storages energy systems like BESS and FESS along with conventional sources like 
Diesel Engine and Fuel Cell systems are integrated into the concept of the hybrid energy system that is 
considered as the microgrid. Many types of research interest in the modeling and control of these kinds of 
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hybrid systems as in [4-7] which have been accomplished to address the control of microgrid using various 
cases of configurations. 

In this paper, a microgrid system is considered to overcome the challenges facing the use of 
renewable generation. This proposed microgrid includes two renewable generation system such as 
photovoltaic (PV) panels and wind generator with diesel engine, and fuel cell system as additional 
generators, the Battery and Flywheel systems are integrated into the microgrid system. From many studies, 
this microgrid configuration is considered as the most important solution for producing electricity in isolated 
areas. However, the major challenge that facing the exploiting the microgrid is to regulate the power and 
frequency fluctuation in an islanded and connected mode of microgrid because of the difficulty of design 
control[8-9]. Autonomous operation of the microgrid, through controlling the system frequency, will enhance 
the performance of the microgrid. 

Due to its simplicity and facility to implement the PID controller becomes the most controller used 
to maintain the system frequency under fluctuations in generation power and load demands[10-11]. To 
achieve the optimal control of frequency in the microgrid system, different optimization techniques are being 
utilized in literature as in [12-17]. These researches applied various optimization techniques like the Genetic 
Algorithm, Mine Blast Algorithm (MBA), Particle Swarm Optimization (PSO), and quasi-oppositional 
harmony search to optimize the controller parameters. Recently few researches have been achieved to present 
frequency control using the Krill Herd technique [18]. However, in this work, the Krill Herd algorithm is 
used to optimizing the PID parameters to control the frequency and power deviations in the proposed 
microgrid. The best-obtained parameters of the PID controller by KH are compared with the others presented 
in [7] using GA. In general, KH is a robust search and powerful optimization technique, able to solving 
global functional optimization problems [19-20]. This method is focused to simulate the behavior of krill 
swarms [21]. The rest of the paper is orderly as bellow. Section 2 describes the proposed configuration of the 
hybrid energy system and its different components models. In Section 3 brief introduction to the PID 
controller and the objective function using in this study. Section 4 explains a review of the Krill Herd 
technique using for the optimization of PID controller gains. Section 5 displays the results analysis and 
comparisons with the performances of the controller using two different optimization techniques and also 
their robustness contra the disconnection of energy storage devices. The paper ended with a conclusion 
in Section 7. 


2. PROPOSED MICROGRID SYSTEM 

The proposed hybrid system model presented by the transfer function is presented in Figure 1. The 
PV array and Wind Turbine Generator (WTG) are considered the principal sources to outfit the load demand. 
The Diesel Engine Generator (DEG) and Fuel Cell, are used as complementary generators to ensure the 
system operation continuity which affected by the interruption nature of renewable sources. The Battery and 
Flywheel devices are added to the system stability. The photovoltaic generation system and the wind energy 
system are combined with other sources and energy storage devices to obtain a more constant power profile. 
This configuration is used in [7]. The microgrid’s parameters are illustrated in Table.1. The typical model of 
the microgrid system utilized in this study is as fellow: [7-10, 22-25] 
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Figure 1. Block of the proposed hybrid energy system 
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Table1l. Proposed microgrid’s parameters 








Component Gain (K) The time constant (s) 
WTG Kwro =1 Twro = 1.5 
7 Tın = 0.04 
EV Koy =I Tic = 0.004 
FC Krc = 0.01 Tre = 4 
DEG Kpgg =0.003 Tore = 2 
BESSS Kpggss =1 Tggss = 0.1 
FESS Kress =1 Tress = 0.1 
— _Kwrc_ _ APwrc 
Gwre = 14+STwrg  4Pwre (D 
_ _Kpv_ _ APpv 
Gpy = 1+STpy AG 2) 
— Krc _ APrc 
Grc(s) = Iste iu (3) 
_ _Kpec_ _ APpgG 
Gpgg(s) = he ae (4) 
—  KpBess _ APsess 
Gsgess (S) = rens A (5) 
_ _ Kress _ APress 
Gress(S) = ina au (6) 


Where Au is the signal control action of the PID controller in feedback to minimize the frequency 
deviation Af. 

The power generation from renewable sources and power demand is modeled in considering the 
small stochastic fluctuations and large deterministic drift [6-7]. 


P= (on. /B(1—G(s)) + B)/B)T =x. (7) 


Where is the power component, P represents the output of wind or solar system and load, P gives 
the value of the power, n is renormalization constant of the generated or load. (x) to correspond the per unit 
(p.u.) level, I is a time-dependent converting signal with a gain that transcribes the sudden change in the 
value for the power output [7, 23]. 

Figure 2 shows the stochastic model of generation powers (Pw, Psol), the total power (Pt), and also 
the demand (PI). Significant fluctuations that can influence the frequency deviation can be remarked. The 
stochastic natures are independent of the control strategy. 
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Figure 2. Realization of the solar, wind and load powers 
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3. PID CONTROLLER STRUCTURE AND OBJECTIVE FUNCTION FOR OPTIMIZATION 

PID controller is used as a specific regulator in loop feedback that is used widely in the industrial 
regulation system. Its standard structure has a ‘‘three-term” controller, which can be modeled using the 
transfer function in its ideal form by (8) or in its parallel form by (9) [11]. 


G(s) = Kp (1 ++ Tps) (8) 
IS 
G(s) = Kp +% 4+ Kps (9) 


Where Kp is the gain of proportionality T;is the time constant of integral, Tp the time constant the derivative, 
K; is the integral gain and Kp is the derivative gain. 

The ‘‘three-term” operations are presented as follows [3, 12]. 

--The proportional term furnishing a global regulation action proportional to the error signal meanwhile the 
all-pass gain. 

-The integral term decreasing steady-state errors over low-frequency using the integral action. 

-The derivative term enhances transient response by high-frequency compensator using differential action. 
Proportional, Integral, and Derivative terms consist of the PID controller actions. 





Plant/process 






Figure 3. Block of PID controller model 


The objective function (J) is considered the optimization problem whole the variables are PID 
controller parameters. It is an integral of the sum square frequency deviation (Af) and the deviation of the 
signal control deviation (Au) as presented by (10) [6-7]. 


Jopt = Jp [WAN + ED Aw ?]dt (10) 


Where w transcribes the relative value of the objectives function (i.e., Integral of Squared Error— 
ISE and Integral of squared Deviation of Control Output—ISDCO), and it is equal to 0.7. Kn = 10* is the 
normalization of scale in ISE and ISDCO. 


4. A REVIEW ON KRILL HERD OPTIMIZATION ALGORITHM 

A new meta-heuristic Algorithm called Krill Herd algorithm is investigated by Gandomi and Alavi 
(2012) that is inspired using the simulation of krill swarm behavior [18] This method is known as robust 
optimization algorithm mimics the behavior of Krill swarms in Krill Herd using for solving difficult 
optimization issues [19]. The krill motion mostly influences the objective function with the distances of each 
krill swarm between food and herd density. The Krill individual’s time position is determined by the 
bellowing three actions [20,25]: 


1. The movement provoked by krill swarms; 
2. Foraging activity; and 
3. Sudden propagation, therefore. 


The provoked movement refers to the density conservation of the herd by each swarm. The 
mathematical formula is as follows: [18-20] 


nee = N” q; + w,N°4 (11) 
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Where N™* is the maximum speed, and æ; is defined as: 
a= qlocal 4 py target (12) 


Wn is the inertia weight of the provoked movement in the interval [0, 1],N°!4 is the last movement provoked, 
q0cal is the local influence given by the neighbors and a'"9° is the target direction effect given by the best 


krill. 


The second movement is influenced by the food present location and the last experience.it can be presented 
as fellow [20]: 


F, = Vbi twp F” (13) 
Where 


Bi = prone Fp pest 


and V; is the speed foraging, wis the inertia weight, F,°'“ is the last foraging movement. The third 
movement is defined as an arbitrary process that is presented with a directional factor and a propagation 
speed. It is expressed by the bellowing formula [20]: 


D; = D™*6§ (14) 


Where D™®* is the maximum diffusion speed and 6 is the random directional vector, and its values are varied 
in [-1, 1]. 


Using various effective parameters of the motion according to the time, the position of a krill 
individual over the interval t to At is given by the bellow formula [20]: 


ax; 


Xit + At) = X,(t) + At— 


(15) 

It may notice that At is one of the most essential parameters and must be correctly determined under 
the problem to be solved. Indeed, this parameter influences as a scale factor for the speed vector. 
Different algorithms inspired by krill can be determined using the movement characteristics of krill 

individuals. The KH algorithm can be summarized by various steps as below [20]. 

1- Data forms: determine the simple limits, resoluteness of Krill Herd algorithm parameter(s). 

2- Initialization: arbitrary initialize the krill swarm into the search environment. 

3- Evaluation of objective function: calculation of each krill swarm objective function with the krill position. 
- Motion calculation: -Movement provoked by other krill swarms, -Foraging movement —Physical 
propagation. 

4- Implement the genetic operators. 

5- Updating: update the krill swarms position in the search environment. 


6- Recurring: return to step 3 until the iteration number is attainted. 
7- End 


5. SIMULATION AND RESULTS 

The proposed system is performed and analyzed under MATLAB/Simulink software. Time-Domain 
analysis of the proposed hybrid system is investigated using KH and GA based PID controller. The Krill 
Herd algorithm is applied to optimize the PID controller’s parameters in the proposed hybrid power system 
illustrated in Figure 1 for frequency control and the obtained results are compared with those obtained by GA 
reported in [7]. Results based on the objective function J solved using the proposed KH is reported. Figure.4 
displays the convergence of the objective function using KH. The convergence characteristics of the Krill 
Herd Algorithm for the PID controller are illustrated in Figure.4. 
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Figure 4. Fitness convergence of the proposed KH algorithm 


Figure 5 shows the frequency and power variations with the signal control deviation. The frequency 
and power deviations are affected by renewable power changing according to the weather conditions. This 
effect is mitigated by the use of the PID controller based Krill Herd algorithm in comparison with the GA. 
The PID controller eliminates the mismatches and enhances also the performance of the system. However, it 
can be shown that the oscillations band is less with KH than with GA. This is more important to facilitate the 
design of the control signal which activates the feedback components such as the FESS, BESS, 
and DEG, etc. 
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Figure 5. Frequency and power variations with the control signal 
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Figure 6 displays the responses of the various generation devices of the proposed microgrid system 
such as BESS, FESS, FC, and DEG under nominal conditions of operation. It can be shown that there are 
fewer fluctuations with PID controller tuning By KH than with the PID controller optimized using GA. 
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Figure 6. Responses of different components of the microgrid system using the best PID-KH and PID-GA 


To limit the frequency variation by the elimination of the unbalance in supply and demand under the 
stochastic change of both the generated power and load demand, the power deviation is controlled by a PID 
controller. The parameters of the controller are optimized using the Krill Herd algorithm and compared with 
the result obtained using GA in [7]. It has been shown that the frequency response using Krill Herd 
algorithm-based PID outperforms the PID controller optimized using the GA algorithm. PID controller tuned 
by the Krill Herd is validated as the better controller over the GA-PID controller. 
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The simulation of the proposed system without battery is presented in Figure 7 and Figure 8. It can 
be observed that the BESS has a high impact on the fluctuation in the frequency and power of the hybrid 
system. This fluctuation caused by the stochastic and intermittent form of wind and solar powers. Also, the 
impact is observed in various power generation by different sources like FC, DEG, and FESS. The results 
show the robustness of the PID controller based KH that minimizes the power and frequency fluctuations in 
the absence of the battery energy storage device. The generated power from different components, as 
reported in Figure 9, shows the robustness of the control scheme against the energy storage disconnecting. 






































T T T T 
gar F Best run values 4 
L Average run values | J 
0.165 - a 
Ay, (0.167 4 
EA L 
3 
x 0.155 4 
ià F 4 
0.15 F 4 
0.145 F m 
fi fi I 1 i 
0 20 40 60 80 100 120 


No. of Iterations 


Figure 7. Objective function convergence of KH algorithm 
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Figure 8. Frequency and power deviations with control signal with and without BESS using PID-KH 
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Figure 9. Power of various components of the autonomous system with without and BESS 


CONCLUSION 


The paper presented an application of a PID controller scheme to regulate the frequency deviation in 


hybrid power system consists of WTG, PV, DEG, FC, BESS, and FESS as illustrated in figure.1 which are 
considered the most promising and sustainable configuration used in the kind of hybrid energy system . Solar 
and wind cause frequency and power oscillations caused by the stochastic nature of this renewable 
generation. Therefore, frequency control is provided by the integration of the PID controller based Krill 
Herd. A comparison with the results offered by GA as reported in [7] is done. These comparisons show 
evidence that the KH outperformed the GA. However, performances of the KH based PID controller is better 
than PID based GA restricting the frequency and power oscillations. Also, the robustness against 
disconnecting the Battery energy storage is tested. 
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